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Isoionic isotope exchange data of 45 Ca 2+ and 32 P0 4 3 ~ in a saturated solution with hydroxylapatite at 25 °C, 
collected by Avnimeleeh, have been reinterpreted on the basis that the processes consisted of three pools: 
hydration shell, surface layer and "recrystallization of crystallites." These processes are regarded as strictly 
separable. This theoretically allows: (i) a quantitative evaluation of the constants involved in the processes; (ii) an 
estimation of the exchange capacities of hydration and surface layers (fast pools) which may be pH-dependent; and 
(iii) a determination of the isoionic exchange rate constants. The constants are essentially the same for Ca 2+ and 
PO-i' 1- at a given pH if the cross-sectional areas of the ions are taken as 23 A 2 and 33 A 2 respectively; they, 
perhaps, increase linearly with the hydrogen ion concentration of the solutions. 

The effects of three consecutive abrupt dilutions of the radioactive ions without changing the concentrations 
(or pH) of the nonradioactive components of the solution are predictable by a mathematical model based on the 
complete reversibility of the two fast pools and the kinetics of the irreversible third pool. 

Key words: Dilution effect, at same pH and concentration; hydration shell, role of; hydroxylpatite, at different 
pHs; isotopic exchange, isoionic; recrystallization, role of; surface layer, role of. 



1 . Introduction 

L 

The determination of the isoionic isotope exchange rates of 
45 Ca 2+ and 32 P0 4 3- in a saturated solution with hydroxylap- 
atite ("bone mineral", OHAp) may have some fundamental 
role in understanding the growth, dissolution, adhesion 
mechanisms and calcium and phosphate homeostasis of bone 
and tooth tissues. An adequate experimental study with 
proper theoretical understanding of the exchange mechanism 
should also throw some light on the surface chemistry of 
hydroxylapatite which is of immense importance in the fields 
of mineral processing and agriculture in addition to dental 
and medical research. 

Exchange experiments have been carried out by a number 
of investigators [1-6]. 1 In some cases, the studies were 
conducted under conditions where other interactions such as 
dissolution, precipitation and hydrolysis were also present 
and, therefore, in such cases the application of exchange 
kinetics is evidently not valid. We have chosen to reinterpret 
the data collected by Avnimeleeh [4-7]. These experiments 
were performed under the conditions of saturation at various 
hydrogen ion concentrations at 25 °C using a homogeneous 



1 Figures in brackets refer to literature references at the end of this paper. 
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nonlabile solid phase which was stabilized by prolonged 
boiling. We have treated Avnimelech's full data, although 
partial treatments have been reported [4, 7]. 

Various theoretical models have been proposed to explain 
the kinetics of the isoionic exchange process with hydroxyl- 
apatite. In general, the isoionic exchange is assumed to 
comprise three [1, 4, 91 J or more [3, 8] processes, each 
characterized by definite kinetic and capacitative properties. 
All models contained two pools, hydration shell and surface 
layer. The third irreversible pool, "recrystallization process," 
may, however, not be unique and may itself be composed of 
two or more pools depending upon the exposed crystallo- 
graphic faces and their energies, etc. Loebenstein's treatment 
[7] is unique. He interpreted the exchange process on the 
basis of Langmuir kinetics for the adsorption of a component 
from solution by a reversible and an irreversible process, 
simultaneously. Loebenstein's approach is valid, but an 
explicit correlation between the assumed adsorption and the 
actual exchange is best approximated only during the early 
stages. 

Treatment of the isoionic exchange by a three-pool proc- 
ess, which we have adopted here, was first considered by 
Neuman et al. [1] The exchange times for the experiments 
considered here were generally less than three hours; it is 
possible that four or more pool processes might have to be 
invoked for longer time intervals. Our approach is different 
from others in that we consider the pool processes as 
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theoretically separable. This allows us to determine quanti- 
tatively the capacities of all three processes and the kinetics 
of the two. The exchange rate constants and the capacities of 
the various pools depend on the hydrogen ion concentrations 
of the solutions. The correlation is, however, not very good 
and it may be attributed to variable contamination of the 
surface layers of hydroxylapatite resulting in its incongruent 
dissolution [10, 11] where Ca to P ratio in solution is not 
5/3. 

Based on our mathematical model, we have been able to 
predict the effect of three consecutive dilutions of the 
radioactive ions. In the dilution process, the concentration or 
the pH of the nonradioactive component is not changed. The 
predictions are generally very good for the first dilutions. For 
the second and third dilutions, where the errors become 
cumulative, the predictions are fair and within the limits of 
experimental deviations of up to 10 and 20 percent respec- 
tively. 

2. Experimental 

The details of the preparation [10] and the experimental 
procedures [4] have been described by Avnimelech. The 
hydroxylapatite was prepared by titrating a calcium hydrox- 
ide slurry at boiling temperature with phosphoric acid in a 
CO 2-free environment. The resultant precipitate was stabi- 
lized by prolonged boiling to minimize recrystallization 
during the subsequent equilibrations and characterized by 
X-ray, infrared, and chemical analyses. Its BET(N 2 ) surface 
area was 26.4 m 2 /g. 



3. Results 

Seven sets of exchange experiments at 25 °C (table 1), 
differing from one another in slurry density of hydroxylapa- 
tite, initial concentration of radioactive tracer ions (i ), 
hydrogen ion concentration (and, therefore, in Ca 2+ and 
PO 4 3_ concentrations), volume of solution, times and fre- 



quency of sampling and time and amount of dilution, are 
examined. The value of solubility product (K s ) determined 
for various runs varies a little. This variation may, however, 
be well within the limits of the incongruent dissolution [10, 
11] of hydroxylapatite. The exchange data for the runs IE10, 
IE8 and IE 13 reproduced here have been reported by 
Avnimelech [4] and quoted by Loebenstein [7]. The exchange 
kinetics, including the dilution effect, for two other runs 
determined by Avnimelech, IE 14 and IE 15, are presented in 
figure 1 for 45 Ca 2+ and in figure 2 for 32 P0 4 3 ~. 

The exchange curves for all sets have three common 
features: (i) a precipitous initial fall within one or two 
minutes, followed by (ii) a gradual decline for up to half an 
hour, and culminating in (iii) a very slow decrease in activity 
for the observational period of up to four hours. For the 
dilution effect, the initial drop in activity in the solution is 
vertical, followed by a slow rise and eventually ending in a 



plateau. 



v! 



4. Discussion 

It seems reasonable [1, 2] to associate the features of the 
exchange curves with the physical processes that may be 
occurring in the system. The initial precipitous drop is 
usually attributed to extremely rapid diffusion between the 
bulk solution and hydration shell; the intermediate decline 
in activity is thought to be linked to the surface exchange of 
ion; [12, 13] and, finally, the very slow decrease may be 
connected with the slow incorporation of ions into the body of 
crystallites which may be primarily due to the recrystalliza- 
tion processes [1, 3, 9]. To facilitate an analysis of the 
exchange kinetics, we consider the pool processes as theo- 
retically separable since the first process is almost instanta- 
neous and the third one, after an hour or so, becomes 
independent of the first two. The kinetics of the second 
process can be determined from the parameters of the first 
and third processes. 



TABLE 1. Solution properties of hydroxylapatite slurries at 25 °C. 



Expt. 


pH 


[Ca 2+ ] 
mol LT 1 X 10 4 


[P0 4 3 -] 
mol L" 1 X 10 4 


(Ca. 2+ f$0 4 3 -f 
(OH-) x 10 59 


Slurry 

density 

g/L 


i b cpm/mL 


designation 


^Ca 


32p 


IE 10 


6.940 


0.667 


0.734 


7.39 


7.212 


25611 


10194 


IE11 


6.718 


0.795 


0.954 


1.64 


3.334 


34425 


17726 


IE15 


6.409 


1.487 


2.144 


3.26 


4.388 


26899 


15011 


IE9 


6.185 


1.921 


3.763 


1.56 


7.352 


18996 


12694 


IE14 


5.817 


3.498 


6.214 


0.27 


7.065 


56415 


39992 


IE8 


5.766 


4.190 


7.980 


0.50 


6.276 


25480 


12610 


IE13 


5.640 


7.469 


14.39 


3.13 


8.193 


56212 


41783 



a The activity products were used to calculate this constant [10]. 
h i is the initial concentration of the radioactive ion. 
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Time, minutes 

FIGURE 1. Exchange kinetics of 4b Ca 2+ with Hydroxylapatite at 25 °C. 
A and B refer to experimental runs IE 14 and IK 15 respectively. 1, 2, and 3 refer to the first, second and third 
dilutions. 
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FIGURE 2. Exchange kinetics o/ 32 P0 4 3 ~ with Hydroxylapatite at 25 °C. 
A and B refer to experimental runs IK 14 and IK 15 respectively. 1, 2, and 3 refer to the first, second and third 
dilutions. 
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The overall radioisotopic exchange mechanism may be radioactive ions, t, present in the solution and, therefore, 
represented as: 

di 

— = k cr Di, (1) 



k h R 

Satd. Soln. ^ Hydration ^ 



dt 



where k cr is the rate constant for the recrystallization process. 



c ^s SL d ^ it 4.- Since S remains constant, eq (1) can be integrated as follows: 

surface — > Kecrystalhzation > n \ / & 



Layer 



where k h , R and k cr are rate constants. The first process is 
almost instantaneous and, therefore, k h is immeasurably 
small and cannot be determined. The two other constants, R 
and A; cr can be determined by our treatment. 

Recrystallization and intracrystalline diffusion. The total 
concentration of calcium and phosphate ions remains con- 
stant in a saturated solution at equilibrium. The dissolution 
of smaller crystals and the growth of bigger ones may, 
therefore, irreversibly remove radioactive ions from the 
solution. Strictly speaking, the mechanism of this slow 
depletion of radioactive ions is primarily due to recrystalli- 
zation and not due to intracrystalline diffusion, although the 
effect may appear the same. The extent of the recrystalliza- 
tion process should be rather small since the rate of loss of 
radioactivity during this period is small and there usually is 
no discernible change in the surface area of hydroxylapatite. 

This is not an equilibrium process and, therefore, the 
McKay equation [14] cannot be applied. It seems reasonable 
to assume that the rate of loss of radioactivity is directly 
proportional to the area of the exposed surface, 5, or the 
slurry density, D, of hydroxylapatite and the concentration of 



In i = —k rr Dt + const. 



(2) 



The constant may be evaluated from the initial condition: 
when t — o, i = i cr ; and it is equal to In i cr . Equation (2) 
now becomes: 



In i = — k rr Dt + In i 



l Kff Lft 

log — = h log 

L 2.303 



(3) 



where i is the initial concentration of a radioactive ion. 
Therefore, the plot of log i/i vs t should be linear and the 
intercept and slope should be equal to log i cr /i () and k cr D/ 
2.303 respectively. These constants for various experimental 
runs, except for IE9 and IE 11 where enough data is not 
available for the recrystallization period, are given in table 2 
together with other statistical parameters. In general, the 
linearity of eq (3) is fairly good and, therefore, the assump- 
tions on which the equation is based should be reasonable. 



Table 2. Recrystallization constants a of hydroxylapatite at 25° C. 



Expt. 
designation 


Slope b 

(k cr D/2, 3) x 10 4 

(min) -1 


Intercept 1 ' 

2 + aog; cr /; ) 


5 C X 10 3 


/yd 


k cr e x 10 4 

L/g. min 


icr 

cpm/mL 


IE10 Ca 


6.92 ± 2.84 


0.875 ± 0.013 


3.90 


4 


2.21 


1922 


IE10 P 


8.28 ± 0.25 


1.054 ± 0.001 


0.29 


4 


2.65 


1154 


IE11 Ca 


(3.07) 


(1.112) 


- 


- 


(2.12) 


(4454) f 


IE11 P 


(3.07) 


(1.289) 


- 


- 


(2.12) 


(3450) 


IE15 Ca 


3.16 ±0.19 


1.279 ± 0.003 


2.64 


8 


1.66 


5116 


IE15 P 


3.05 ± 0.26 


1.589 ± 0.005 


3.78 


8 


1.61 


5831 


IE9 Ca 


(4.66) 


(1.252) 


- 


- 


(1.46) 


(3394) 


IE9 P 


(4.66) 


(1.575) 


- 


- 


(1.46) 


(4753) 


IE14 Ca 


4.11 ± 0.68 


0.970 ± 0.005 


3.29 


4 


1.34 


5266 


IE 14 P 


3.70 ± 0.74 


1.254 ±0.007 


3.33 


4 


1.20 


7175 


IE8 Ca 


2.51 ±0.10 


1.434 ± 0.001 


1.89 


10 


0.92 


6917 


IE8 P 


1.89 ± 0.23 


1.682 ± 0.003 


5.02 


14 


0.69 


6065 


IE13 Ca 


2.49 ± 0.88 


1.047 ±0.007 


6.24 


5 


0.69 


6264 


IE 13 P 


2.74 ± 0.53 


1.349 ± 0.005 


4.64 


7 


0.76 


9326 



a Based on eq (3) of the text. 

b Least-squares estimate of the coefficient ± its standard error. 

c Standard error of regression. 

d Number of measurements. 

e The slurry densities instead of the related surface areas of hydroxylapatite are used to normalize the values of k c 

f Values in parenthesis are obtained from the linear dependence of A: rr on pH and best estimates of i cr values. 
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FIGURE 3. Variation of the exchange rate constant for recrystallization process with pH of the solution at 
25 °C: A refer to 45 Ca 2+ exchange and O refer to 32 P() 4 3 ~ exchange. 
The straight line is obtained by linear regression. 



The values of k rr and i cr may be obtained from the slopes 
and intercepts respectively (table 2, columns 6 and 7). The 
slurry densities instead of the surface areas were used to 
normalize the values of k cr . 

It is observed that k cr , both for calcium and phosphate 
ions, is linearly related to the pH of solution: 



k cr = a x + 6 t pH. 



(4) 



The plot is shown in figure 3. Linear regression gives the 
values of constants: a x = -(6.18 ± 0.82) X 1CT 4 , b 1 = 
(1.23 ± 0.13) X 10~ 4 and 5 = 2.08 X 10~ 5 , where s is the 
standard error of regression. An analysis of variance indi- 
cates [15] that, for a 95 percent confidence level, no gain is 
made in assigning separate linear equations to Ca 2+ and 
PO4 3- . For IE9 and IE1 1, the values of k cr are obtained from 
the constants of eq (4) and the i cr values are the best 
estimates. 

Surface Layer Exchange. This exchange may be viewed as 
occurring between the hydration shell, H, and the surface 
layer, L, of hydroxylapatite and may be represented as: 

HX* + LX ^± LX* + HX, 



where X * is a radioisotope exchanging with a normal isotope 
X. Since we have assumed that the processes are separable, 
it is evident that the radioisotopic concentrations may be 
taken as i e at the end of the first process or at the beginning 
of the second one and i cr at the end of the second process or 
at the beginning of the third one. Now, applying the McKay 
equation, [14, 16] we may easily obtain that: 



log(i' - i cr ) =log(i e - i cr ) 



R(a + b)t 
2.303ab 



(5) 



where a (mol/L) is the total concentration of exchanging ion 
in solution and b (mol/L) is the concentration of ion on the 
surface of hydroxylapatite which is in equilibrium with one 
litre of the solution. The term b is a composite term and may 
be expressed as: b = N L calc D, where /V L calc (mol/g) is the 
calculated concentration of total exchangeable ion on the 
surface (V L calc = S/Na, see footnotes to table 3) and/) (g/L) 
is the slurry density of hydroxylapatite. 

The plot of log(i — i cr ) vs t should be linear and the 
intercept and slope should be equal to log(i e — i cr ) and R(a 
+ b)/2.303ab respectively. These constants for various 
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Table 3. Constants associated with surface layer exchange a at 25 °C. 



Expt" 
designation 


Slope 
R(a + b) 10 2 

2.303a/3 


Intercept 

log {i e - i cr ) 


s c X 10 2 


N 


R x W 5 
mol/L. min 


i e 

cpm/mL 


R* x 10 5 
mol/L. min 


IE 10 Ca 


8.60 ± 0.61 


3.288 ± 0.033 


3.41 


6 


1.26 


3863 


1.15 


IE 10 P 


7.19 ± 0.40 


3.026 ± 0.035 


5.00 


7 


1.13 


2216 


0.98 


IE11 Ca 


9.22 ± 0.25 


3.613 ± 0.012 


1.45 


7 


1.50 


8575 


1.33 


IE11 P 


7.85 ± 0.33 


3.300 ± 0.022 


3.42 


9 


1.42 


5445 


1.13 


IE15 Ca 


6.19 ± 0.51 


3.500 ± 0.073 


8.76 


5 


1.80 


8281 


1.42 


IE15 P 


4.02 ± 0.13 


3.311 ± 0.025 


3.55 


6 


1.45 


7879 


0.80 


IE9 Ca 


8.19 ± 0.19 


3.290 ± 0.021 


1.66 


3 


3.19 


5344 


2.43 


IE9 P 


5.80 ± 0.15 


3.030 ± 0.018 


1.40 


3 


3.63 


5825 


1.66 


IE 14 Ca 


5.49 ± 0.25 


3.482 ± 0.033 


4.63 


6 


3.51 


8301 


3.52 


IE 14 P 


6.42 ± 0.67 


3.356 ± 0.063 


5.56 


4 


5.52 


9443 


5.26 


IE8 Ca 


4.91 ± 0.55 


3.471 ± 0.064 


8.72 


7 


3.51 


9878 


2.49 


IE8 P 


4.89 ± 0.39 


3.171 ± 0.069 


9.54 


8 


4.58 


7548 


2.60 


IE 13 Ca 


8.08 ± 0.49 


3.305 ± 0.044 


6.08 


7 


9.40 


8281 


9.54 


IE13 P 


5.06 ± 0.14 


3.161 ± 0.016 


2.02 


7 


7.21 


10774 


6.30 



a Based on eq 5 of the text where b = /V L cak D — SD/Ncr and N £ calc (mol/g) is the calculated concentration of total exchangeable ion on the surface; S 
is the surface area (26.4 m 2 /g) and D is the slurry density (g/L, Table I) of hydroxylapatite; N is the Avogadro's number and cr is the effective area of the 
exchangeable ion (23 A 2 for Ca 2+ and 33 A 2 for P0 4 3_ ). 

b All points obtained before two minutes have been neglected. 

c Standard error of regression. 
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Figure 4. Log (i - i cr ) vs lime for 45 Co 2+ at 25 X. 


A and B refer to experimental runs IE10 and IE1 1 respectively. The straight lines are obtained by linear regression. 
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experimental runs are given in table 3 together with other 
statistical parameters and figures 4 and 5 show the plots for 
two runs. In general, the linearity of eq (5) is fairly good. 
The experimental points do not, generally, fall on the straight 
line (Hgs. 4 and 5) before two minutes have elapsed. This is, 
possibly, because it takes at least that much mixing time for 
the solution to reach uniformity after a radioactive ion is 
injected into it. Our analysis shows that the exchange 
between surface layer and hydration shell or bulk solution is 
complete within half an hour after which the linearity of eq 5 
is not maintained. Similar equilibrium times were observed 
for nickelous [12, 13], cupric [12, 13] and zinc [17] ions 
exchanging with hydroxylapatite surface. 

Radioisotopic exchange experiments at the point of zero 
charge [5] (pH = 8.5) and exchange experiments with other 
ions [12, 13, 17] with hydroxylapatite yield effective areas of 
23 A 2 and 33 A 2 for calcium and phosphate ions respectively. 
These values, together with the BET(N 2 ) surface area of the 
hydroxylapatite yield the calculated concentration of a total 
exchangeable ion ([V L calc )on the surface. 



The constants, R and i e , can be derived from the slope 
and intercept of eq (5), respectively (table 3, columns 6 and 
7). It is observed that /?, both for calcium and phosphate 
ions, is linearly related to the hydrogen ion concentration of 
solution: 



R = a 2 + b 2 [H + \ 



(6) 



The plot is shown in figure 6. Linear regression gives the 
values of constants :a 2 = (0.82 ± 0.46) X 10" 5 , b 2 = (0.274 
± 0.037) X 10~ 5 and s = 1.09 X 10" 5 wh 



wnere s is 



the 

standard error of regression. An analysis of variance indi- 
cates [15] that, for a 95 percent confidence level, no gain is 
made in assigning separate linear equations to Ca 2+ and 
P0 4 3 ~. The expected identity of/? values obtained from the 
exchange of either calcium or phosphate ions also demon- 
strates that the assumed ratio [5] (23/33) of the cross- 
sectional areas of the two ions is reasonable. 

The dependence of/? on [H] + in eq (5) does not necessarily 
imply that the intercepts should also be so dependent. In 




Time, minutes 

Figure 5. Log (i - i cr ) vs time for 32 P0 4 3 ~. 

A and B refer to experimental runs IE10 and IE11 respectively. The straight lines are obtained by linear regression. 
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fact, a casual perusal of the intercept values (table 3, column 
3) shows that they are fairly constant 3.30 ± 0. 15. 

Capacities of Hydration Shell and Surface Layer. The 
capacity of the hydration shell (N H ) can be easily derived as: 



[Bulk Solution] [Hydration Layer] 



^v 



M*< 



M* H 



[Surface Layer] 
M* 



(10) 



N„ = 



a i n - l 



D 



(7) 



where symbols have their usual designations. The total 
experimental capacity of hydration shell and surface layer is: 



where the numerators respectively designate the concentra- 
tions of total exchangeable ions in solution, in hydration 
layer and on hydroxylapatite surface, and the denominators 
give the respective concentrations of radioactive ion (If*) in 
various pools. In the present case eq (10) may be represented 



(N H + N L ) = 



a I i n — ir 



D 



(8) 



The experimental capacity of the surface layer ( / V L ) is the 
difference of eqs (8) and (7) 



N,=- 



I) 



( l lgr \ _ / l () 
i cr ] \ 






(9) 



DN, 



This can also be obtained from the following basic condition 
of equilibrium of a radioisotope in various reversible pools: 



\i i cr ) i c1 



T- (ID 



It is evident that eq (11) is identical to eq (9). 




*. 



10 12 14 16 18 20 22 24 



[H^xlO 7 , mol/L 



FIGURE 6. Variation of exchange rate constants for surface layer with [H^J of the solution at 25 °C: #, ■ refer to Ca 2+ and O, □ refer to PO4 
The straight lines are obtained by linear regression. 
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The values of the capacities of various pools for different 
runs are given in table 4. The various capacities do not 
change very much with hydrogen ion concentration up to 
pH = 6; the change is continuous and dramatic after that. 
One such plot is shown in figure 7. It may be that the 
variation {N (H or L or H+D vs [H + ]} is parabolic. The data is, 



however, scarce and scattered, and other possibilities cannot 
be ignored. The scattering of data may be related to the 
variable contamination of the surface layers of hydroxylapa- 
tite [10, 11]. Probably, the surface concentrations of constit- 
uent ions are best obtained at the point of zero charge of the 
substrate [5]. 



TABLE 4. Capacities of hydration shell and surface layer at 25 °C. 





Ca 2+ capacity 3 , mol/g 


PO4 3 capacity 11 , mol/g 


(IS, +N H ) 


Expt. 
designation 


[Ca 2+ ][£] 

S 
x io 4 


[Ca 2+ ][,4 - B] 

S 
x 10 4 


x 10 4 


[poriB] 

s 

x 10 4 


[P0 4 3 -][/t - B] 

S 
x 10 4 


(N L + A/„) x 
10 4 


for Ca 2 + 
divided by 

Wl + Nh) 
for PO4 3 - 


IE10 


0.521 


0.619 


1.140 


0.366 


0.431 


0.797 


1.430 


IE11 


0.719 


0.886 


1.605 


0.645 


0.539 


1.184 


1.356 


IE15 


0.762 


0.681 


1.443 


0.442 


0.327 


0.769 


1.876 


IE9 


0.668 


0.534 


1.201 


0.604 


0.252 


0.855 


1.405 


IE14 


2.870 


1.939 


4.809 


2.845 


1.178 


4.023 


1.195 


IE8 


1.054 


0.738 


1.792 


0.853 


0.519 


1.372 


1.306 


IE13 


5.277 


1.993 


7.269 


5.055 


1.057 


6.112 


1.189 


Mean ± standan 


deviation 












1.39 ± 0.23 



a A = (i — i cr )/i cr and B = (i — i^li f ,. N H and N L are the capacities of the hydration shell and surface layer respectively. 
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Figure 7. Variation of total capacity of reversible pools ( N H + N L ) with [H + ] of the solution at 25 °C: • refers to Ca 2+ and O refers to PO4 
The curves are best parabolic fits 18 . 
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The mean of the ratios of the total capacities (hydration 
shell 4- surface layer) of calcium ions to phosphate ions is 

1.39 (table 4, last column). The theoretical ratio is 33/23 = 
1.43. The two values agree within the limits of experimental 
error. 

By substituting the experimental value of the total capacity 
of surface layer, N L (table 4, columns 3 and 6), for N L calc in 
the slope of eq (5) (table 3, column 2), the true value of the 
surface exchange rate constant (R *) can be calculated (table 
3, column 8). The values of R * are, however, not much 
different from those of/?. The rate constant/?*, also, varies 
linearly with [H + ] (fig. 6). Linear regression yields the values 
of constants: intercept = (0.35 ± 0.59) X 10~ 5 , slope = 
(0.260 ± 0.047) X 10" 5 and standard error of regression = 

1.40 X 1(T 5 . 

Dilution Effect . The effect of dilution on the concentration 
of radioactive ions, without changing the concentrations (or 
the pH) of nonradioactive components of the solution, may 
be divided into reversible and irreversible parts: 



i'), 



(12) 



where i lt is equilibrium concentration of a radioactive ion, 
corrected for the irreversibly removed ions, after the first 
dilution at a given time, i r l is the equilibrium concentration 
for reversibly exchanging ions after the first dilution, i cr has 
already been defined, and i t is the concentration of radioac- 
tive ion present in the solution after a given recrystallization 
time (eq (3), table 2). 



The constant i r l can be obtained from the following 
equilibrium condition: 



[Solution] / 
[Solution] j f 



[Surface Layer] 
[Surface Layer] , 



[Hydration Layer] x l 
[Hydration Layer] , f 



(13) 



where (,*) refers to the beginning and (j f ) to the end of the 
first dilution. The concentration of a radioactive ion after the 
dilution in any one of the pools may be easily calculated 
from the appropriate material balance equation. Considera- 
tion of any one equality between the three terms in eq (13) is 
sufficient to derive i t \ We shall derive it from the equality 
between the first and second terms, although it may be easier 
to derive from the equality between the first and third terms. 
Considering the equality between the first and second terms, 
eq (13) may be physically represented as: 



M {1 ° 



*)] 



»-H(^] 



(14) 



where V ^ is volume of the solution before the first dilution 
and V] f after the dilution (table 5, column 2). A simple 



Table 5. Radioactive exchange and dilution effect at 25 °C. 







First dilution 


Second dilution 












[i r l - (i cr -i % )]% 


Third dilution 


Expt. 


w 


cpm/mL 












designation 


vf 


Expt 


Calcd. 




i r ",c 


om/mL 




i r ra cpm/ 
mL 




Expt. 


Calcd. 


Expt. 


Calcd. 


IE 10 Ca 


505/705 


1685 


1623 


645/945 


1540 


1559 


890/1290 


1370 


1476 


IE10 P 


" 


950 


933 


" 


850 


881 


" 


720 


816 


IE11 Ca 


425/625 


4075 


4088 


571/871 


3575 


3744 


817/1217 


3200 


3345 


IE11 P 


" 


3050 


3072 


" 


2390 


2709 


" 


1950 


2320 


IE15 Ca 


490/690 


4140 


4090 


640/940 


3760 


3717 


890/1290 


3350 


3294 


IE15 P 


" 


4310 


4194 


" 


3480 


3541 


" 


2800 


2904 


IE9 Ca 


534.5/834.5 


2750 


2789 


798.5/1398.5 


2375 


2379 


— 


— 


_ 


IE9 P 


" 


3625 


3636 


" 


2850 


2721 


- 


— 


_ 


IE14 Ca 


500/700 


4350 


4419 


650/950 


3950 


4206 


900/1300 


3300 


3938 


IE 14 P 


" 


5400 


5885 


" 


4250 


5374 


" 


3300 


4789 


IE8 Ca 


662/1062 


5120 


5071 


— 


— 


— 


— 


— 


_ 


IE8 P 


" 


4225 


4090 


— 


— 


— 


— 


— 


_ 


IE13 Ca 


490/690 


5200 


5502 


645/945 


4550 


5170 


900/1300 


3675 


4765 


IE13 P 


" 


7000 


7744 


" 


5350 


6900 


" 


4050 


5990 



1 The recrystallization times (in minutes) in order of experimental runs are: 85, 35, 221, 84, 141, 245 and 143. 
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manipulation of eq (14) yields: 



L + L 



C4 



(15) 



The same result is also obtained if no hydration layer is 
assumed. The various terms of eq (15) are known and, 
therefore, i r l can be calculated. 

The amount of radioactive ion (i*), present in the solution 
after a given recrystallization time, may be calculated from 
eq (3) (slopes and intercepts are given in table 2). Therefore, 
the contribution of the irreversible component, (i cr — i*), can 
be determined. 

The corrected values ( i lt ) of the equilibrium concentrations 
of an ion after the first dilutions at any given time can now be 
calculated from eq (12) and are presented in table 5 (column 
4). The experimental values are the best estimates from the 
enlarged graphs (i versus t) at a given time (generally, the 
time when maxima are attained and this is the time for which 
i f is calculated). The agreement between the experimental 
and calculated values of i lt is very good, considering the 
scattering of the experimental data on one hand and the 
limitations of theoretical treatment on the other. 

It is difficult to apply a correction (due to the irreversible 
removal of ions) to the equilibrium concentration of an ion 
after the second or third dilution, since it is not possible to 
identify the experimental quantities with theoretical terms. 
And, therefore, no attempt is made to apply the correction 
term for second and third dilutions. 

The concentration of a radioactive ion after the second 
dilution, considering the exchange as completely reversible, 
may be best calculated from the following equality: 



[Solution] 2 
[Solution] 2 



[Hydration + Surface Layer] 2 l 
[Hydration + Surface Layer] 2 f 



(16) 



The concentrations in the second term can, again, be 
calculated from the material balance equation and eq (16) 
may be represented as: 



v 2 f" V J 



(17) 



constituting eq (14). A simplification of eq (17) gives: 






(18) 



Similarly, the concentration of a radioactive ion after the 
third dilution, considering the exchange as completely re- 
versible, may be calculated as: 



(19) 






where the new symbols have designations analogous to those 



The values of i r I] and j r m are given in table 5. The agreement 
between the calculated and experimental values of i r n and 
i r !I1 appears reasonable for most runs. If the error due to the 
irreversible removal of radioactive ion from the solutions 
could be estimated, a better agreement probably would 
result. 

A better correlation for the exchange experiments, con- 
ducted at various hydrogen ion concentrations, might result 
if the ineongruent dissolution of hydroxylapatite due to the 
variability in the composition of the surface layers were 
overcome by rejecting the initially equilibrated solutions two 
or three times before starting the experimental runs. 



6. Glossary 

total concentration (mol/L) of exchangeable 
ion, 

total concentration (mol/L) of exchangeable 
ion on surface of hydroxylapatite which is 
in equilibrium with one liter of solution, 
slurry density (g/L) of hydroxylapatite, 
concentration (cpm/mL) of radioactive ion 
in solution at time t, 

concentration (cpm/mL) of radioactive ion 
in solution at the end of the second process 
(diffusion from hydration shell to surface 
layer) or at the beginning of the third one 
(recrystallization process), 
concentration (cpm/mL) of radioactive ion 
in solution at the end of the first process 
(diffusion from bulk solution to hydration 
shell) or at the beginning of the second one 
(diffusion from hydration shell to surface 
layer), 

equilibrium concentrations (cpm/mL) for re- 
versibly exchanging radioactive ions after 
first, second and third dilutions respec- 
tively, 



D 



i r \ i r { \ i r m 
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i lt = equilibrium concentration (cpm/mL) of ra- 

dioactive ion, corrected for the irreversibly 
removed ions, after first dilution at time "t" ', 

k cr = rate constant (L/g. min) of exchangeable ion 

for recrystallization process, 

N H — capacity (mol/g) of hydration shell for ex- 

changeable ion, 

N L , N L calc = experimental and calculated capacity (mol/ 
g) of total exchangeable ion on surface, 

/?,/?* = experimental and "true" surface exchange 

rate constant (mol/L. min) of ion, 

S = surface area (m 2 /g) of hydroxy lapatite, 

Vi\ V2 1 , ^V — initial volumes (mL) of solutions before first, 
second and third dilutions respectively, 

Vi *, V2 f , ^3 f — final volumes (mL) of solutions after first, 
second and third dilutions respectively. 

7. References 

[1] Weikel, J. H., Neuman W. F., and Feldman, I., J. Am. Chem. Soc, 

76, 5202 (1954). 
[2] Neuman, W. F., and Neuman, M. W., The Chemical Dynamics of 

Bone Mineral, University of Chicago, Chicago, 111., 1958. 
[3] Pak, C. Y. C, and Bartter, F. C, Biochim. Biophys. Acta, 141, 401 

(1967). 



[4] Avnimelech, Y., Israel J. Chem., 6, 375 (1968). 

[5] Kukura, M., Bell, L. C, Posner, A. M., and Quirk, J. P., J. Phys. 

Chem., 76,900(1972). 
[6] Kukura, M., Bell, L. C, Posner, A. M., and Quirk, J. P., Soil Sci. 

Amer. Proc, 37,364(1973). 
[7] Loebenstein, W. V., J. Coll. & Interface Sci., 36, 234 (1971). 
[8] Zimens, K. Arkiv for Kemi (Mineralogi 0. Giologi), 20A, 1 (1945). 
[9] Edgington, D. N., Radiat. Res., 25, 257 (1965). 
[10] Avnimelech, Y., Moreno, E. C, and Brown, W. E., Solubility and 
surface properties of finely divided hydroxyapatite, J. Res. Nat. 
Bur. Stand. (U.S), 77A, (Phys. and Chem.), No. 1, 149-155 (Jan.- 
Feb. 1973). 
[11] Smith, A. N., Posner, A. M., and Quirk, J. P., J. Colloid Interface 

Sci., 48,442(1974). 
[12] Misra, D. N., Bowen, R. L., and Wallace, B. M., J. Colloid Interface 

Sci., 51,36(1975). 
[13] Misra, D. N., and Bowen, R. L., J. Biomed. Matls. Res., 12, 505 

(1978). 
[14] McKay, H. A. C, Nature, 142, 997 (1938). 

[15] Afifi, A. A., and Azen, S. P., Statistical Analysis: A Computer Oriented 

Approach, (Academic Press, Inc., New York, N.Y., 1972), pp. 

108-113, 146-151. (Thanks are due to H. Argentar who pointed 

out this method and reference to the author. ) 

[16] Frost, A. A., and Pearson, R. G., Kinetics and Mechanism, (John 

Wiley and Sons, Inc., New York, N.Y., 1953), p. 179. 
[17] Misra, D. N., and Bowen, R. L., to be published. 

[18] Mood, A. M., and Graybill, F. A., Introduction to the Theory of 
Statistics, (2nd ed.), (McGraw Hill Book Co., New York, N.Y. 
1963). 



406 



